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ABSTRACT 


An  analytical  formulation  of  the  bed  shear  stress  coef- 
ficient inside  the  surf  zone  is  derived  using  the  concept 
of  radiation  stress.  A truncated  Rayleigh  p.d.f.  is  used 
to  describe  the  wave  field  inside  the  surf  zone  and  provides 
the  input  to  calculate  the  variation  of  wave  energy  and  long- 
shore current  as  a function  of  wave  height,  water  depth  and 
distance  to  shore.  The  wave  set-up  is  approximated  using  a 
sinusoidal  wave  solution.  Field  measurements  of  longshore 
current  and  waves  within  the  surf  zone  are  used  to  calculate 
the  bed  shear  stress  coefficient.  The  data  consist  of  647 
data  points  selected  from  LEO  program  and  62  data  points 
from  Ingle  (1966)  observations,  all  taken  along  the  Southern 
California  coast.  Frequency  distributions  and  statistics 
are  calculated  for  the  bed  shear  stress  coefficient.  A mean 
bed  shear  stress  coefficient  to  two  significant  decimal 
places  is  found  to  be  0.01. 
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INTRODUCTION 


A.  STATEMENT  OF  THE  PROBLEM 

It  is  known  that  when  sea  waves  or  swell  approach  a 
straight  coastline  at  an  oblique  angle  a mean  current  is 
generated  parallel  to  the  shoreline,  see  Figure  I.  Such 
longshore  currents  are  of  prime  importance  for  both 
coastal  engineering  and  for  aiding  in  the  strategic  planning 
of  Naval  inshore  warfare  operations. 

An  accepted  theory  of  longshore  currents  on  plane 
beaches  is  developed  in  terms  of  the  momentum  flux  due  to 
the  waves  directed  down  coast  being  balanced  by  the  shear 
stress  associated  with  the  mean  flow.  The  formulation  of 
the  bed  shear  stress  requires  the  specification  of  a bed 
shear  stress  coefficient.  The  purpose  of  this  thesis  is 
the  determination  of  the  bed  shear  stress  coefficient  to  be 
used  in  the  longshore  current  formulas. 

The  study  will  also  help  in  the  analysis  of  sediment 
transport.  The  shear  stress  does  work  on  the  bottom  in 
moving  sediments.  Several  authors  have  formulated  sediment 
transport  in  terms  of  the  bed  shear  stress  which  in  turn 
requires  an  appropriate  bed  shear  stress  coefficient. 

B.  HISTORICAL  REVIEW 

Inman  and  Quinn  (1952),  using  the  momentum  approach  for 
the  prediction  of  longshore  current  by  Putnam,  Munk  and 


> CURRENT 


Figure  I 


Definition  of  Longshore  Current  Variables 


Taylor  (1949),  showed  that  in  order  to  fit  theory  with  obser- 
vations, the  bed  shear  stress  coefficient  must  be  permitted 
to  vary  with  the  longshore  velocity  over  a wide  range  of 
3 1/2  orders  of  magnitude. 

Bretsc hne i de r (1954)  found  that  the  spectral  limitations 
of  wave  growth,  under  the  action  of  steady  wind  in  shallow 
water  with  a typical  sandy  bottom,  suggested  a value  for  the 
friction  coefficient  of  between  0.01  and  0.02.  Also,  he 
found  that  the  observed  damping  of  swell  propagating  over  a 
smooth,  leveled,  impermeable  sea  bed  was  consistent  with  a 
value  of  the  coefficient  of  between  0.034  and  0.097. 

Longuet-H i gg  i n s (1970),  using  the  concept  of  radiation 
stress,  developed  a relationship  for  prediction  of  the 
theoretical  maximum  longshore  current  just  inside  the  break- 
ing and  proposed  a friction  coefficient  of  the  order  of  0.01. 
He  concluded,  on  the  basis  of  the  finding  of  Bretschne i der 
(1954),  Prandtl  (1952)  and  Nikuradse’s  experiment  with 
roughened  pipes,  that  is  was  not  "...unreasonable  to  expect 
a friction  coefficient  of  the  order  of  0.01." 

Table  I was  taken  from  Sonu  (1975);  it  summarizes  some 
values  of  the  friction  coefficient  proposed  by  various  inves- 
tigators. The  values  reported  were  obtained  from  measure- 
ments outside  the  surf  zone  or  from  laboratory  experiments. 

It  can  be  seen  from  this  table  that  the  range  of  values  is 
relatively  wide  and  the  test  conditions  varied. 
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TABLE  I 


Bottom  Friction  Coefficients 
Proposed  by  Various  Investigators 


Friction 

Coe  fficient 

Wave 
Height 
(meters ) 

Wave 

Period 

(sec) 

Test 

Con d i t ion  s 

Authors 

0.01 

Arbitrary 

Arbitrary 

Sha 1 low  water 
steady  state 
wave  generation 

Bretschnei der 
( 1954, a) 

0.030-0.089 

0.23-0.51 

2.88-3.96 

Gulf  of  Mexico; 
depths  3. 4-5.2  m 
slope  0.00035- 
0.00-41 

Bretschnei der 
( 1954, b) 

0.030-0.040 

2 

10 

Niigata,  Japan; 
depths  2.25- 
2.75  m 
slope  0.018 

Kishi  (1954) 

0.01 

5.4* 

8.4* 

Osci 1 lating  water 
channel,  turbu- 
lent boundary 
layer 

Jonsson  (1966) 

0.01  -0.40 

0.002- 

0.100 

0.88-2.58 

Wave  flume, 
laminar  boundary 
layer 

Iwagaki  and 
Tsuchiya( 1966) 

0.03  -0.18 

1 .77-2.47 

9.1-15.5 

Hiyshizu,  Japan; 
depths  13-10  m 
slope  0.0060 

Iwagaki  and 
Kakinuma( 1 966) 

0.03  -0.15 

1.05-1.60 

7.4-12.5 

Takahama,  Japan; 
depths  10-7  m 
slope  0.0057 

Iwagaki  and 
Kakinuma( 1 966) 

0.09  -0.50 

(in72irn)-4~  20 

U:  nearbottom  velocity 
n:  ripple  height 

Wave  flume  study; 
derived  from 
energy  dissipa- 
tion in  sand 
ripple  vortices 

Tunsta 1 1 and 

1 nman  ( 1975) 

Note:  *Equi valent  values  at  10  m depth 


C.  OBJECTIVES  OF  THE  STUDY 


The  bed  shear  stress  coefficients  previously  determined 
are  based  on  a very  limited  set  of  field  data  or  on  labora- 
tory studies  which  used  as  a model  simple  sinusoidal  waves 
which  are  not  typical  of  the  randomness  found  in  nature. 

The  objective  of  this  study  is  to  analyze  existing  sets  of 
field  observations  obtained  in  the  surf  zone  and  by  using 
the  best  available  theory  attempt  to  determine  a reasonable 
value  of  the  bed  shear  stress  coefficient.  For  this  purpose 
a fairly  large  data  set  obtained  for  the  Channel  Island 
Littoral  Environment  Observations  (LEO)  Program  was  used  as 
weli  as  a set  of  observations  by  Ingle  (1966)  taken  at  various 
locations  along  the  Southern  California  coast.  It  is 
expected  that  the  data  obtained  and  the  theory  applied  will 
ultimately  contribute  to  the  establishment  of  a reasonable 
value  of  the  bed  shear  stress  coefficient  and  to  a more 
accurate  prediction  of  the  longshore  current  velocity  across 
the  surf  zone. 
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THEORY 


I I . 

A.  INTRODUCTION 

Several  models  have  been  proposed  for  the  distribution 
of  the  longshore  current  velocity  across  the  surf  zone  on  a 
plane  sloping  beach.  The  so  I ut i on , wh i ch  uses  pure  sinusoids 
to  describe  the  waves  and  no  lateral  shear  stress,  gives  a 
velocity  distribution  which  is  triangular  shaped  with  both 
a peak  velocity  and  a discontinuity  at  the  breaker  point 
E C Bowen,  I 969 ) , Thornton  (1969)  and  Longuet-H i gg i ns  (1970)3. 
This  is  unreasonable  since  there  are  no  discontinuities  in 
nature.  A second  model  including  lateral  shear  stress  tends 
to  smooth  out  the  discontinuity  at  the  breaker  line  and  pro- 
duces a smoother  velocity  distribution  with  the  maximum 
velocity  occurring  closer  to  shore.  However,  no  criterion 
to  predict  an  optimum  lateral  shear  stress  coefficient  is 
as  yet  available.  This  introduces  an  added  complexity  to 
the  problem. 

A random-sea  model  developed  by  Collins  (1972)  circum- 
vents the  difficulty  of  the  lateral  shear  stress  coefficients 
and  allows  the  statistical  input  of  the  sea  state  as  de- 
scribed by  a Rayleigh  distribution.  Figure  2 compares  the 
velocity  distribution  resulting  from  the  various  models:  the 
non  lateral  stress  model,  the  iateral  shear  stress  model  for 
a coefficient  equal  to  0.4  ( Longuet-H i gg  i ns , 1970),  and  the 
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random-sea  model.  The  velocities  are  referenced  to  the  peak 
velocity,  Vm,  of  the  non  lateral  shear  stress  model,  which 
according  to  Longuet-H i gg  i ns  (1970)  is  given  by 


Vm 


5-nq 

T6C^ 


s i na 


and  V is  the  longshore  component  of  the  mean  current  velocity 
Bowen  (1969),  Thornton  (1969)  and  Long uet-H i g g i n s (1970) 
attributed  the  generation  of  longshore  currents,  due  to  an 
oblique  wave  approach,  to  the  longshore  component  of  the 
momentum  flux  (radiation  stress)  of  the  water  waves.  The 
calculation  of  the  wave-induced  longshore  current  velocities 
and  changes  in  mean  water  level  requires  the  specification 
of  the  radiation  stresses  as  a function  of  the  location  and 
wave  properties  in  the  nearshore  region. 

B.  WAVE  SET-UP  INSIDE  THE  SURF  ZONE 

As  waves  approach  the  coast  and  shoal  there  is  a change 
in  the  momentum  flux  of  the  waves  which  is  balanced  by  a 
change  in  the  mean  water  level.  Outside  the  surf  zone  there 
is  a set-down  while  inside  the  surf  zone,  after  breaking, 
there  is  a set-up  or  superelevation  of  the  water  level.  The 
wave  set-up  is  important  because  both  the  local  wave  height 
and  speed  are  functions  of  the  total  local  water  depth  which 
is  unknown.  The  change  in  mean  water  level  required  to 
balance  the  excess  momentum  flux  of  the  waves  must  be 
determined  first.  A convenient  form  of  the  x-component 
(shoreward)  of  the  momentum  flux  equation  integrated  over 
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depth  and  averaged  in  time,  derived  by  Longuet-Higgins  and 
Stewart  (1962),  for  describing  the  wave  set-up  is  given  by 


+ 35  * 0 


( I ) 


which  says  that  the  change  of  excess  momentum  flux  due  to 
wave  action  ("radiation  stress")  is  balanced  by  a change  in 
the  mean  water  level.  It  is  assumed  in  the  derivation  that 
the  net  local  mass  flux  perpendicular  to  shore  is  zero  so 
that  there  is  no  contribution  from  the  mean  motion  to  the 
momentum  flux  perpendicular  to  the  shore,  and  that  the  mean 
stresses  are  negligible. 

Inside  the  surf  zone  it  is  assumed  that  the  radiation 
stress  tensor  can  be  expressed  in  terms  of  the  energy  and 
wave  speed  in  the  same  form  as  in  shallow  water.  This  implies 
that  the  breaking  waves  are  of  the  spilling  type  and  that 
even  under  breaking  waves  the  water  particle  motion  retains 
much  of  its  organized  character  as  described  by  linear  wave 
theory.  Using  the  shallow  water  approx imat i on. that  the  group 
velocity  equals  the  phase  velocity  and  that  the  angle  cf 
incidence  equals  zero,  the  radiation  stress  term  reduces  to 


3 .2-.  .-<2 


Sxx  = 76  p9Y  (h+n)' 


(2) 


Substituting  (2)  into  (I)  and  integrating  gives  the  mean 
water  elevation  of  the  form 


n = K(h.-h)  + n.  , 
b b 


(3) 
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where 


i 


I 


i 


K 


I + 


_l 

1 9 

378Y 


(4) 


and  the  mean  water  elevation  at  breaking 


is  given  by 


(5) 


is  a single  breaker  height  which  corresponds  to  the 
significant  wave  height  observed  at  the  breaking  position. 
For  this  study,  it  was  assumed  that  the  waves  followed 
the  break i ng  i ndex. 


Y 


0.78  , 


derived  from  a modified  solitary  wave  theory  by  Munk  (1949). 
The  total  local  depth  of  water,  D,  is  obtained  by  combining 
the  local  mean  water  super  elevation  and  the  local  depth, 

D = n + h . (6) 


A sinusoidal  description  of  the  waves  was  used  to  solve 
for  the  wave  set-up  in  order  to  get  a closed  form  analytical 
solution  and  to  circumvent  the  difficult  numerical  solution 
of  equation  (I)  required  by  the  random-sea  model.  Collins 
(1972)  compared  the  sinusoidal  solution  to  the  random  sea 
solution  as  shown  in  Figure  3.  The  effects  of  the  random-sea 
model  is  to  smooth  the  waves  set-up  curve.  The  magnitude  of 
the  difference  between  the  two  solutions  is  very  small  but 
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wave  sot-up  on  offshore  slopes  of  1:25  and  1:50,  broken 
(.•suits  for  pet  iodic  waves  having  the  same  energy  content. 

Figure  5.  Comparison  of  Wave  Set-up  Solutions. 


the  percent  difference  can  be  large  as  the  depth  of  water 
approaches  zero.  The  total  depth  of  water  is  important  in 
prescribing  the  breaker  point  or  the  I imits  of  integration 
on  the  Rayleigh  distribution  for  the  random-sea  model. 

In  the  calculation  of  longshore  currents  using  the  random 
sea  model  it  is  the  area  under  the  Rayleigh  distribution 
that  is  used  so  that  small  errors  on  the  limits  generally 
cause  only  even  sma I ler  errors  in  the  area.  Hence,  it  is 
felt  using  sinusoidal  wave  descriptions  to  calculate  wave 
set-up  is  a reasonable  approximation. 

C.  WAVE  FIELD  INSIDE  THE  SURF  ZONE 

A description  of  the  wave  field  is  required  in  the  long- 
shore current  calculation  because  knowledge  of  it  is  needed 
for  specifying  the  horizontal  water  particle  velocities  and 
for  determining  the  longshore  component  of  the  radiation 
stresses  in  an  irregular  wave  field.  Inside  the  surf  zone 
the  waves  are  unstable  and  the  fluid  motion  loses  some  of 
its  ordered  character;  but  Thornton  (1976)  points  out  that 
most  of  the  water  particle  motion  in  the  body  of  the  fluid 
is  coherent  with  the  surface  and  can  be  considered  wave-in- 
duced and  not  turbulent,  particularly  for  spilling  type 
brea  kers . 

In  this  study,  a truncated  Rayleigh  distribution  as  shown 
in  Figure  4 is  used  to  give  a statistical  description  of  the 
wave  field  as  described  by  Collins  (1972)  and  Battjes  (1974). 
The  basic  assumption  is  that  at  each  depth  a limiting  breaker 


height  can  be  defined  which  cannot  be  exceeded  by  the  indi- 
vidual waves  of  the  random  field,  and  that  those  wave 
heights  which  in  the  absence  of  breaking  would  exceed  the 
breaker  height  are  reduced  by  breaking  to  the  value  of  the 
local  breaker  height.  That  is,  the  energy  corresponding  to 
the  height  in  excess  of  the  local  breaker  height  is  assumed 
to  be  dissipated.  The  limiting  breaker  height  decreases  as 
depth  decreases. 

In  describing  the  Rayleigh  distribution,  a fictitious, 

or  reference,  local  energy  per  unit  area,  denoted  as  E r , is 

defined.  The  reference  energy  density  refers  to  that  energy 

density  that  would  exist  if  breaking  had  not  occurred  nad 

accounts  for  shoaling  and  refraction  transformation.  Battjes 

(1974)  also  defines  reference  wave  heights  Hr  and  their  mean 
2 

square  value  H which  is  related  to  Er  according  to 

Er  = F pg  ^7"  * (7) 

The  reference  wave  heights  are  assumed  to  be  Rayleigh  distri- 
buted . 

The  Rayleigh  wave  height  distribution  is 

clipped  at  H = H s in  accordance  with  the  assumption  that  the 
height  of  a breaking  wave  equals  the  local  breaker  height, 

Hs,  in  order  to  obtain  an  approximation  to  the  actual  wave 
height  distribution.  Then,  the  mean  energy  per  unit  area  at 
a fixed  point,  taking  account  of  breaking,  is  calculated  from 

E = JL  pg  H2  . (8) 
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The  variance  is  calculated  from  the  pd*  of  H, 


where 


p ( H ) dH 


p ( H ) = j d[l-exp(-H2/Hr2)] 
o 

+ H 2 exp ( -H  2/H  2 ) , 
s s r 

H2  - J"  H2d[l-exp(-H2/H^T)] 

o 

+ Hs2  exp(-Hs2/Hr2)  , 


H2  = [l-exp(-H  2/H2)]  H2  . 

r s r r 


(9) 


(10) 


(II) 


(12) 


The  clipped  Rayleigh  distribution  implies  that  all  waves 
from  H s to  infinity  that  were  previously  larger  than  Hs  now 
are  reduced  to  the  same  height,  H . Therefore,  the  total 
probability  (percent)  of  waves  having  the  height  Hs  is  given 
by 

OO 

/ p ( H ) dH  . 
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The  contribution  to  the  variance  is  given  by 


H 


s 


00 


p ( H ) dH  = H 2 
s 


exp  ( - 


9 


s 

which  is  the  term  on  the  right  of  equation  (II). 
local  breaker  height,  which  inside  the  surf  zone 
to  be  given  by 


H i s the 
s 

i s assumed 


Hs  = YD  . (13) 

The  local  mean  wave  height  H can  be  expressed  in  a similar 
manner  in  terms  of  Hs  and  Hr  by  means  of  the  clipped  Rayleigh 
distribution, 

/H 

HdCl-exp(-H2/hr7)]  + H s exp(-Hs2/iT7) 
o (14) 

where  again  the  term  on  the  right  represents  the  percent  of 
waves  greater  than  Hs  in  the  original  distribution  which  now 
have  the  height  H . Integrating  (14)  gives 

/T 

H = Hr  erf (Hs/Hr  ) (15) 

rms  ' rms 


in  which 


rms 


and  the  error  function  being  defined  as 


erf ( p ) 


— f 

/r  J 


exp(-t  )dt  . 


The  error  function  was  calculated  using  the  rational  approxi 
mation  of  Abramowitz  and  Stegun  (1965), 


erf(p)  = I — C ( a . t + a„t2  + a^t"5)  exp(-p2)3  + £(p) 


where 

t 

z 

P 


I / ( I +zp  ) , 

.47047, 

H /(H  2) 1 /2, 
s r ' 


a j = .3480242, 
a 2 = .0958798, 
a3  = .7478556. 


The  largest  error  using  this  approximation  is 
e ( p ) <_  2.5  x I0"3  . 

In  the  observations  used  for  comparison  with  the  theory, 
the  breaker  height  is  measured  visually.  It  is  assumed 
that  an  observer  visually  measures  the  significant  breaking 
wave  height  defined  as  the  average  of  the  highest  one-third 
fraction  of  the  wave  heights.  The  difficulty  in  applying 
this  definition  to  the  present  problem  is  that  the  defini- 
tion applies  to  a point  measurement  or  a statistically  homo- 
geneous (spatial)  wave  field  and  in  this  problem  the  waves 
are  defined  as  varying  spatially  as  they  shoal  shoreward. 
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In  order  to  define  the  significant  wave  height  for  a 
spatially  varying  ( non homogeneou s ) wave  field,  it  is 
assumed  that  the  observer  measures  waves  when  spatially 
one-third  of  the  waves  have  broken;  hence,  the  reference 
wave  height  can  be  specified  from  the  clipped  Rayleigh  dis- 
tribution 

exp  ( -H . 2/H  ^ ) = 0.333 
b r 

and  H 2 = - H 2/ Ln ( 0 . 333 ) . (16) 

r b 


D.  LONGSHORE  CURRENT  VELOCITY 

The  derivation  of  the  longshore  current  starts  with  the 
y-momentum  flux  equation 


3S 

+ T -T,  - 0 

ax  b z 


(17) 


Where  the  lateral  shear  stress,  x^,  is  neglected  (17)  reduces 
to 


as 


3x 


*1  + xb  = 0 , 


(18) 


which  says  that  the  cnange  of  y-component  (longshore) 
momentum  flux  due  to  waves  in  the  x-direction  is  balanced  by 
the  bottom  stress,  x^,  in  the  y-direction.  Assuming  that 
the  amplitude  of  the  wave  motion  |Uw|  is  much  greater  than 
the  mean  current  velocity,  V,  then  (Thornton,  1969) 


xb  = p Cf  | Uw | V . 


(19) 


r 
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The  excess  of  momentum  flux  of  the  waves,  or  "radiation 
stress"  component  is  given  by 


S = E s i na  cosa  , 
xy  C 

which  inside  the  surf  zone  reduces  to 

S = E s i na  cosa  ( 20 ) 

xy 

under  the  assumption  that  = C for  shallow  water.  Combin- 
ing equation  (3)  and  (12)  gives  the  mean  energy  per  unit 
area  as 

E = -g-  pg  Hr2[l-exp(-Hs2/Hr2)3  (21) 
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from  equations  (3),  (5),  (6),  and  (13)  Hs  can  be  expressed 


as 


H = H.  (K  - -fr-)  + h ( x ) (y-YK)  . 

S D 10 


(23) 


Application  of  Snell's  law  allows  the  local  breaker 
angle,  a,  to  be  expressed  in  terms  of  the  known  breaking 
angle  afa  and  the  breaking  celerity  Cb, 


sina  = 77—  sin  a,  . 
Cb  b 


(24) 


Using  the  shallow  water  approximation  for  wave  speed 


C = (gh ) 


1/2 


and 


H. 


r - fn  \l/2 
Cb  ' (g  ~ > 


then 


sina  = ( h r^-)  * ^ sina. 

Hb  b 


(25) 


Hence,  from  equations  (23)  and  (25)  it  is  seen  that  both 


H and  a are  now  expressed  as  functions  of  h,  which  in  turn 


is  a function  of  x. 

The  bed  shear  stress  coefficient  is  determined  by  com- 
bining equation  (18)  and  (19), 

3S 


n x y 

cf  = 7x  p | Uw | v * 


(26) 
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The  mean  horizontal  water  particle  velocity  amplitude  is 
expressed  using  linear  theory  (Battjes,  1974), 


Uw 


■ £ £ 
it  D ’ 


(27) 


where  H is  given  by  equation  (15), 

The  change  in  the  radiation  stress  is  given  by 


3S 

9x 


r i/o  Yh  sin  a, 

21  , l(JL)1/2  si™, t(h. 


2 7 7 

~ 2yhsin  a. 


2 H 


S(  I - 


) + 


+ smct4cosa  pgS(Y-YK)Hs  exp(-Hs2/Hr2)  , (28) 


where 

E = i pg  ir7[l-exp(-Hs2/(77)] 


and 


. y .1/2 

si  na  cosa  = (rr~) 

Hb 


2 . 2 
yh  sin  a 


s i na,  ( h - 
b 


b)  I / 2 
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DATA 


A.  LEO  DATA 

Data  from  the  Littoral  Environment  Observation  (LEO)  pro- 
gram established  by  the  Coastal  Engineering  Research  Center 
was  used  in  this  study.  In  the  LEO  program,  nearly  simul- 
taneous observations  of  breaker  conditions  (height,  period, 
angle  of  approach  and  type),  local  winds,  longshore  currents, 
foreshore  slope,  width  of  the  surf  zone  and  rip  currents  were 
made  daily  during  the  period  under  consideration.  The  long- 
shore current  was  determined  by  observing  the  direction  and 
measuring  the  distance  parallel  to  shore  that  a dye  packet 
injected  into  the  surf  zone  traveled  in  one  minute.  Appendix 
A provides  the  set  of  instructions  followed  during  the  obser- 
vations. 

The  data  used  for  this  study  cover  a period  from  May  1972 
to  September  1975  and  refer  to  stations:  5703,  5706,  5707, 
5713,  5714  and  5715,  located  within  the  confines  of  Point 
Mugu  Naval  Air  Station,  60  miles  northwest  of  Los  Angeles, 
California  (location  6 in  Figure  5). 

From  these  sta t i ons, 4 , 632  observations  were  considered 
of  which  only  647  data  points  were  used  in  the  analysis.  The 
following  criteria  were  discussed  to  eliminate  observations 
which  were  not  consistent  with  the  application  of  the  theory 
or  were  simply  erroneous: 
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I . Rip  Currents 


It  was  beyond  the  scope  of  this  study  to  account  for 
any  modifying  effect  of  the  longshore  current  system  by  rip 
currents.  Hence,  all  observations  noting  the  presence  of 
a rip  current  were  systematically  deleted. 

2 . Angle  of  Wave  Approach 

The  theory  used  in  the  derivation  of  the  formulas 
employed  in  this  study  assumes  that  the  angle  between  the 
direction  of  wave  approach  and  the  depth  contour  must  be 
different  from  90  degrees  in  order  for  a longshore  current 
to  be  generated.  Hence,  all  observations  in  which  the  wave 
direction  was  reported  as  being  perpendicular  to  the  shore- 
line were  neglected. 

3 . Wind 

Shepard  and  Inman  (1950)  suggested  the  importance 
of  the  wind  in  generating  longshore  current;  they  also  indi- 
cate that  it  is  difficult  to  separate  the  wind  generated 
current  contribution  from  the  current  generated  by  the  waves. 
Thus,  observations  where  the  wind  speed  was  reported  as 
being  greater  than  ten  miles  per  hours  were  not  considered. 

4 . Foreshore  S I ope 

Observations  where  the  foreshore  slope  was  reported 
as  being  greater  than  ten  degrees  were  neglected  since  such 
large  values  are  not  consistent  with  what  is  usually  observed 
on  the  beaches  under  consideration. 

5 . Wave  Period 

Arbitrarily,  to  keep  the  study  restricted  to  sea  and 
swell  of  relatively  short  period,  all  observations  where  the 
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wave  period  was  reported  as  greater  than  20  seconds  were 
neg I ected . 

6 . Doubtf u I Data 

All  observations  in  which  the  reported  data  were 
considered  to  be  incorrect  due  to  either  mistakes  of  the 
observer  or  the  typist,  such  as  longshore  currents  in  ex- 
cess of  six  feet  per  second,  direction  of  approach  greater 
than  180  degrees,  distance  of  dye  injection  greater  than 
600  feet,  etc.,  were  systematically  rejected. 

B.  SOURCES  OF  ERROR 

Considering  the  interest  in  longshore  currents,  it  is 
somewhat  surprising  that  there  are  relatively  few  sets  of 
adequate  field  measurements  of  longshore  currents  and  the 
simultaneous  wave  parameters  in  the  surf  zone.  After  a 
search  of  the  literature  it  was  concluded  that  little  has 
been  achieved  for  devising  electronic  equipment  designed 
for  gathering  longshore  current  and  associated  wave  informa- 
tion on  a routine  basis.  Hence,  as  in  the  case  of  the  LEO 
data,  most  of  the  observations  must  rely  on  the  good 
judgment  and  personal  abilities  of  the  observers.  This 
introduces  a subjectivity  factor  which  ultimately  affects 
the  final  results. 

I . Brea  ker  Angle 

Galvin  and  Nelson  (1967)  suggested  that  the  variable 
most  difficult  to  measure  with  necessary  accuracy  is  the 
angle  of  wave  approach  or  wave  direction.  Galvin  and  Savage 
(1966)  suggested  that  when  using  a visual  compass  referenced 
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to  a baseline  to  measure  the  breaker  angle,  the  errors  may 
easily  be  +_  two  degrees,  leading  to  a relative  error  which 
is  very  large  for  small  breaker  angles  but  which  decreases 
as  the  breaker  angle  increases. 

In  the  LEO  observations  a protractor  was  used  for 
determining  the  breaker  angle  as  shown  in  Appendix  A.  This 
system  is  completely  visual  using  the  unaided  eye  to  estab- 
lish the  perpendicular  to  the  shore  and  introduces  a greater 
human  factor.  Hence,  it  is  a good  assumption  to  attach  an 
accuracy  less  than  that  suggested  by  Galvin  and  Savage  to 
such  measurements. 

2 . Beach  Slope  and  Surf  Zone  Width 

In  describing  longshore  currents  which  flow  within 
the  surf  zone,  accurate  knowledge  is  required  of  the  beach 
profile  including  both  the  beach  slope  and  the  width  of  the 
surf  zone.  The  approach  used  for  the  LEO  data  was  to  assume 
a plane  beach  inside  the  breaker  line.  The  nearshore  sub- 
aqueous slope  was  computed  using  the  observed  surf  zone 
width  and  the  observed  breaker  height.  Thus  an  uncertainty 
factor  for  the  beach  was  introduced.  Referring  to  LEO  obser- 
vation instructions  in  Appendix  A,  the  observation  of  the 
surf  zone  width  "...  is  based  upon  the  judgment  of  the  obser- 
ver; man-made  or  natural  features  in  the  surf  (e.g.,  a pier) 
may  aid  in  this  observation."  Again,  a subjectivity  factor 
is  involved. 

The  observed  foreshore  slope  could  not  be  used  because 
it  proved  to  be  unrepresentative  of  the  beach  profile  within 
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the  surf  zone 


The  calculated  water  depths  inside  the  break- 


> 
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er  line  using  this  slope  were  systematically  greater  than 
the  calculated  breaker  depths. 

3 . Wave  Period 

Galvin  and  Nelson  (1967)  point  out  that  under  favor- 
able conditions  the  wave  period  can  be  measured  with  reason- 
able consistency  in  the  field  by  visual  observation.  Although 
this  parameter  was  not  used  directly  in  the  computation  it 
is  interesting  to  notice  that  their  suggestion  agrees  quite 
well  with  the  LEO  observations  since  the  range  of  periods 
found  fall  into  the  expected  values  for  the  shoreline  under 
consi de ration. 

4 . Breaker  Height 

In  the  LEO  program,  the  breaker  height  observation 
is  based  solely  on  the  judgment  of  the  observer.  Known 
dimensions  of  natural  or  man-made  features  on  the  shoreline 
or  in  the  surf  zone  are  used  as  references  for  estimating 
the  wave  breaking  height.  Galvin  and  Savage  (1966)  sug- 
gested a relative  error  in  breaker  height  measurement  of 
+_  25  percent.  They  arrived  at  this  figure  by  comparison  of 
breaker  height  measurements  made  with  pressure  gages, 
oscillographs  and  visual  observations,  although  the  measure- 
ments were  not  made  simultaneously.  Hence,  in  the  light 
of  their  finding  it  can  be  concluded  that  at  least  the  same 
error  should  be  expected  in  the  LEO  data  in  which  the  obser- 
vations are  solely  visual. 
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5 . Longshore  Current 


In  the  LEO  program,  the  current  speed  was  determined 
by  using  a dye  as  a tracking  agent.  This  also  adds  an  uncer- 
tainty factor  due  to  the  diffusivity  characteristic  of  the 
dye. 

C.  INGLE  DATA 

A set  of  62  field  observations  made  along  the  Southern 
California  coast  (Fig.  5),  taken  from  Ingle  (1966)  and  per- 
sonal communication,  were  selected  using  the  same  criteria 
used  for  selecting  the  LEO  data.  Despite  the  size  of  Ingle 
sample,  about  ten  percent  the  size  of  the  LEO  data,  its 
analysis  is  important  since  the  Ingle  observations  are  more 
accurate  than  the  LEO  observations.  Thus,  the  Ingle  results 
serve  as  a reference  comparison  to  the  results  obtained  using 
the  LEO  data. 

The  parameters  H^,  period,  ab  and  V were  taken  directly 
from  a summary  appendix  in  Ingle  (1966).  The  beach  profiles 
and  the  distance  shoreward  from  breaker  zone  in  which  the 
longshore  current  velocities  were  recorded  also  were  avail- 
able. However,  the  positions  of  the  breaker  were  not  avail- 
able. Thus,  the  parameter  hb  was  obtained  from  the  relation- 
ship hb  = 1.28  Hb;  and  the  local  depth  below  still  water 
level,  h,  and  the  beach  slope  were  scaled  out  from  the  beach 
profiles  presented  in  the  publication. 

It  should  be  pointed  out  that  in  the  Ingle  data  the  break- 
er heights  were  measured  by  sighting  on  either  a graduated 
pole  held  at  an  approximate  still-water  line  and  the  horizon, 
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a graduated  pole  held  in  the  zone  of  breaking  waves,  or  a 
piece  of  cardboard  with  a slit  and  graph  paper  along  one 
edge.  Breakers  less  than  2 feet  in  height  were  estimated 
while  standing  in  the  breaker  zone.  In  the  LEO  program, 
breaker  height  observations  were  based  solely  on  the  judg- 
ment of  the  observer  on  the  shore.  For  measuring  the 
breaker  angle,  Ingle  observers  used  a Brunton  compass  while 
standing  in  the  surf  zone,  supplemented  by  sights  taken 
from  positions  elevated  above  the  beach;  LEO  observers  used 
a protractor,  as  shown  in  Appendix  A,  with  the  observer  on 
the  beach.  For  measuring  other  parameters,  both  LEO  and 
Ingle  observers  used  essentially  the  same  techniques.  It 
is  important  to  mention,  that  in  the  case  of  the  Ingle 
observations,  the  beach  slope  and  the  position  of  the  obser- 
vations in  the  surf  zone  were  better  than  those  of  the  LEO 
observations  since  in  the  former  an  ordinate  and  abscissa 
arrangement  of  wooden  stakes  al lowed  workers  to  position 
themselves  in  the  surf  zone;  also  most  people  involved  in 
Ingle  observations  were  well  trained  personnel. 
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IV.  RESULTS 


r 


Equations  (3),  (5),  (6),  (12),  (15),  (16),  (21),  (22), 
(23),  (25),  (26),  (27)  and  (28)  presented  in  the  theory  sec- 
tion were  used  with  the  parameters  H^,  a^,  x^,  x and  V from 
LEO  and  Ingle  field  observations  to  solve  for  the  bed  shear 
stress  coefficient.  The  coefficient  calculated  is  based  on 
data  at  specific  locations  within  the  surf  zone  and  not  for 
( mean  conditions.  It  should  be  mentioned  that  since  the  co- 

efficient is  not  determined  by  direct  measurement,  it  there- 
fore  not  only  reflects  bed  shear  stress,  but  also  any  errors 
and  uncertainties  in  measurement. 

The  significance  of  the  assumption  that  the  observer 
measured  waves  correspond  to  a clipped  Rayleigh  distribution 
when  spatially  one-third  of  the  waves  have  broken,  used  for 
computing  the  reference  wave  height  H , was  tested  using 
other  assumptions.  It  might  just  as  logically  be  argued 
that  the  significant  wave  height  might  correspond  to  the 
point  where  half  the  waves  have  broken,  in  which  case 

exp(-Hs^/H^  = 0.5  . 

This  assumption  was  used  for  computing  new  values  for  the 
coefficient.  The  relative  difference  between  the  coeffi- 
cients thus  calculated  and  the  original  ones  was  determined; 
the  variability  was  found  to  be  of  the  order  of  ten  percent. 


which  is  relatively  low.  Hence,  it  can  be  said  that  the 
bed  shear  coefficient  calculations  are  not  very  sensitive 
to  the  assumed  definition  of  the  significant  wave  height. 

A.  LEO  DATA 

Figure  6 depicts  the  frequency  distribution  of  values 
obtained  for  the  coefficient  and  selected  statistics  of  the 
distribution.  The  variability  in  coefficient  values,  as 
represented  by  the  standard  deviation  of  the  distribution, 
is  a measure  of  the  consistency  of  the  calculation  of  the 
coefficient  from  the  field  observations.  The  mean  of  the 
distribution  is  0.008,  while  the  standard  deviation  is  0.010. 
This  suggests  there  is  a large  spreading  in  the  results. 
However,  it  should  be  noticed  that  more  than  90  percent  of 
the  calculated  values  fall  between  0.001  and  0.020  and  that 
the  distribution  has  less  spread  than  a Gaussian  distribu- 
tion for  the  same  standard  deviation. 

B.  INGLE  DATA 

Figure  7 depicts  the  frequency  distribution  of  values 
obtained  for  the  coefficient  and  selected  statistics  of  the 
distribution.  The  mean  of  the  distribution  is  0.014  and 
the  standard  deviation  0.01  I , which  suggests  again  a large 
spreading  in  the  results.  However,  95  percent  of  values 
lie  between  0.001  and  0.030  and  again  the  distribution  has 
less  spread  than  a Gaussian  distribution  for  the  same 
standard  deviation. 
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COEFFICIENT 


Mean 

0.008 

Skewness 

3.694 

Variance 

0.00009 

Ku  rtos i s 

1 7 . 880 

Std.  Deviation 

0.00972 

Minimum 

0.00  1 

Coef . Variation 

1 . 1 9004 

Maximum 

0.080 

Figure  6.  Frequency  Distribution  of  Coefficient  Values 
(LEO  Data). 


COEFFICIENT 


Mean 

0.014 

Skewness 

1 . 584 

Variance 

0.0001 

Ku  rtos i s 

2.649 

Stq.  Deviation 

0.0112 

Minimum 

0.00  1 

Coef.  Variation 

0.8109 

Max  i mum 

0.056 

Figure  7 . 


Frequency  Distribution  of  Coefficient  Values 
(Ingle  Data). 


C.  COMPARISON  OF  RESULTS 


Despite  the  difference  in  sample  size  between  the  LEO 
and  Ingle  data  used  for  the  calculations,  some  comparisons 
can  be  made.  A simple  way  of  comparing  both  sample  results 
is  by  looking  at  their  mean  and  standard  deviation.  The 
mean  and  standard  deviation  corresponding  to  the  values  of 
the  coefficient  for  the  Ingle  data  are  both  larger  than  the 
values  obtained  for  the  LEO  data.  There  is  a relative  dif- 
ference of  75  percent  between  the  mean  of  the  coefficient 
values  of  the  two  samples;  but,  the  relative  difference 
between  the  two  standard  deviations  Is  only  ten  percent. 
This  says  that  the  distribution  values  for  the  coefficient 
in  both  samples  is  nearly  the  same,  although  for  the  Ingle 
data  the  values  for  the  coefficient  were  somewhat  larger 
and  with  more  spread  than  for  LEO  data  which  might  be  ex- 
pected for  the  smaller  sample  size. 

It  was  stated,  when  comparing  both  sets  of  data  the 
Ingle  observations  were  more  accurately  taken  and  more 
reliable  than  the  LEO  observation.  Hence,  the  results  ob- 
tained with  the  Ingle  data  would  be  expected  to  be  better 
than  the  results  obtained  with  the  LEO  data.  To  test  if 
there  is  any  statistical  difference  between  the  two  sets 
of  data,  a hypothesis  test  about  the  two  means  obtained  was 
made.  The  central  limit  theorem  states  that,  if  x is  the 
mean  of  a random  sample  of  size  n taken  from  a population 
having  the  mean  y and  the  finite  variance  a2,  then 

z = hu. 

o/n 
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is  the  value  of  a random  variable  whose  distribution  function 


approaches  that  of  the  standard  normal  distribution  as  n 
The  variances  of  the  population  are  unknown,  but  since  both 
samples  are  fairly  large,  it  is  justifiable  to  approximate 
the  population  variances  with  the  samples  variance.  Thus,  a 
test  statistic  can  be  stated  as 


z 


The  hypothesis  to  be  tested  is  the  null  hypothesis,  y j — y L = 0, 
against  the  alternative  hypothesis  > 0,  where  y repre- 

sents the  mean  of  the  population.  The  evaluation,  for  the 
data  available,  of  the  z statistic  was  found  to  be  equal  to 


z 


0.014  - 0.008 
0.0001  , 0.00009, I / 2 

62  647  ’ 


For  a level  of  significance  of  0.001  the  z statistic  for  the 
normal  distribution  is  3.49.  Since  the  value  obtained  for 
the  test  statistic  is  larger  than  the  critical  value  of  3.49, 
the  null  hypothesis  is  rejected  with  great  confidence;  and  it 
can  be  concluded  that  the  difference  between  both  means  is 
statistically  significant  and  cannot  be  attributed  to  chance. 
Therefore  +he  results  obtained  with  Ingle  data  are  better 
than  the  results  obtained  with  the  LEO  data. 
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D.  CORRELATION  WITH  INDEPENDENT  VARIABLES 


Attempts  were  made  to  correlate  the  calculated  coeffi- 
cients with  the  independent  variables,  breaker  type  and  wave 
period,  which  were  recorded  in  the  field  but  which  were  not 
used  directly  in  the  computations.  Analysis  showed  nothing 
conclusive  regarding  the  correlation  of  the  coefficient  to 
the  breaker  type  since  the  distribution  of  breaker  types 
among  the  data  was  very  uneven;  the  spill/plunge  type  repre- 
sented 72  percent  of  the  data  and  the  spilling  type  20  per- 
cent. Table  II  shows  some  selected  statistics  of  the  dis- 
tribution of  coefficient  values  for  various  breaker  types. 

TABLE  II.  Selected  Statistics  for  Distribution  of 
Coefficient  According  to  Breaker  Type 
(LEO  Data) 


No . 

Observation 

Sp i 1 1 i ng 

1 30 

Plunging 

27 

Surg i ng 
25 

Sp i 1 1 /P 1 unge 
454 

Mean 

0.0090 

0.0075 

0.0098 

0.00788 

Variance 

0.0001 

0.00004 

0.00015 

0.00009 

Std.  Dev. 

0.0107 

0.00609 

0.01204 

0.00946 

Coef.  Var. 

1 . 1 880 

0.0150 

0.23321 

1.  20318 

Range 

0.0790 

0.0230 

0.05300 

0.07500 

Minimum 

0.0010 

0.0010 

0.00100 

0.00100 

Max  i mum 

0.0800 

0.0240 

0.05400 

0.07600 

S kewness 

3.720 

1.3971 

2.90456 

3.74730 

Ku  rtos i s 

17.60 

1 .0553 

7.38989 

18.8567 

A simple  linear  regression  between  the  calculated  coef- 
ficients and  the  observed  period  gave  the  selected  statistics 
of  Table  III. 


TABLE  III.  Correlation  of  Coefficient  with 
Wave  Period  Statistics 


Cor  re  I at i on  ( R ) 

Std.  error  of  estimate 
R squared 
Significance 
I ntercept 
Slope 


LEO  Data  INGLE  Data 


05370 

- .08479 

00970 

0 1120 

00288 

.00719 

08574 

.2  56  1 7 

01088 

.01892 

00023 

- .00042 

The  negative  sign  of  the  correlation  coefficient  indi- 
cates that  there  is  an  inverse  relationship;  that  is,  the 
value  of  the  coefficient  tends  to  become  smaller  as  the 
period  increases.  However,  this  relationship  is  very  weak 
as  indicated  by  the  absolute  value  of  R which  in  both  cases 
is  much  smaller  than  one.  This  result  is  not  surprising 
since  waves  in  shallow  water  become  non-d i spers i ve  or  in- 
variant of  period. 
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V.  CONCLUSIONS 


An  analytical  solution  for  the  bed  shear  stress  coeffi- 
cient was  derived  using  the  concept  of  radiation  stress. 

The  best  theory  for  calculating  the  variation  of  wave 
energy  and  longshore  current,  and  the  resulting  bed  shear 
stress  coefficient,  was  to  use  the  truncated  Rayleigh  p.d.f. 
for  the  statistical  description  of  the  wave  field  inside  the 
surf  zone.  A sinusoidal  approximation  of  the  waves  was  used 
to  calculate  the  wave  set-up.  Calculations  of  the  coeffi- 
cient were  made  by  using  suitable  sets  of  data  obtained 
during  the  LEO  observation  program  and  Ingle  (1966)  observa- 
tions along  the  Southern  California  coast. 

Variability  in  the  results  obtained  for  the  coefficient 
values  were  expected  due  to  subjectivity  and  uncertainties 
in  the  techniques  used  in  the  data  collection.  This  is  the 
first  test  of  the  bed  shear  stress  coefficient  using  fairly 
large  setsof  field  measurements  within  the  surf  zone.  Even 
with  the  uncertainties  involved,  the  analysis  resulted  in  a 
fairly  good  agreement  between  the  mean  of  the  calculated 
coefficients  in  this  work  and  the  values  obtained  by  various 
investigators  for  the  bed  shear  stress  coefficient  for  dif- 
ferent test  conditions  and  outside  the  surf  zone. 

It  was  shown  that  the  dependence  of  the  coefficient  on 
the  wave  period  is  negligible,  in  agreement  with  the  assump- 
tion that  waves  inside  the  surf  zone  are  non-d i spers i ve  or 
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period  invariant.  Since  one  of  the  biggest  differences 
between  Pacific  and  Atlantic  coast  waves  is  the  period,  it 
may  be  concluded  that  the  calculated  coefficient  is  not 
ocean  dependent. 

Since  it  was  initially  concluded  that  Ingle's  data  was 
of  higher  quality  than  the  LEO  data,  it  is  assumed  the 
coefficient  values  using  Ingle's  data  is  therefore  more 
reliable.  In  any  event,  the  mean  value  of  the  two  data 
sets  are  the  same  to  two  significant  decimal  places.  There- 
fore, it  is  concluded  that  a reasonable  value  for  the  bed 
shear  stress  coefficient  within  the  surf  zone  is  0.01. 


APPENDIX  A 

LITTORAL  ENVIRONMENT  OBSERVATIONS 

rw  fens  No.  113-72-8  Mar  72  has  bem  designed  far  keypunching  onto  acnputer  cards  directly  (■■ell  runters  above  each  box  represent  card  ooluei  rasters) . 

It  is  Lsi.uisaiilsil  a p«tcil  be  used.  All  data  should  be  recorded  carefully  and  legibly.  Errors  should  be  corrected  by  first  erasing  erroneous  data  as 
write  overs  usually  produce  illegible  data.  Make  ronarks  as  neaessary  an  t*ie  form  but  record  only  data  in  the  boies  provided.  All  observations  Bust  be 
m«ie  at  the  sane  point  on  the  bsath  every  time  (in  front  of  the  reference  pole) . 

STATION  lOEfflFTCAriCN: 

Cacti  site  m the  "Littoral  Bwirormsnt  Cbeervation"  study  has  been  assigned  a rujneric.il  code  consisting  of  5 digits.  The  first  two  digits  define  the  state 

or  territory  in  which  the  site  is  located  and  the  nrsuunq  3 digits  define  the  particular  bearfi  or  par*  within  the  state  or  territory.  A space  hae  also 

bwi  provided  to  write  in  the  rune  of  the  particular  beach  cr  park  at  wtuch  the  observation  is  taking  place. 

3OT: 

Indicate  in  the  spaces  provided  the  year,  .Tenth  and  day  on  which  an  observation  is  made. 

TOC: 

Iniicate  the  tine  at  which  tto  observations  are  being  made.  The  24-hour  systan  of  recording  tine  has  been  selected  in  order  to  eliminate  any  aon fusion  be- 
tween AM  and  ?M.  The  hour  "00"  refers  to  midnight.  "07"  to  7:00  AM,  "13"  to  1:00  W,  etc. 

SUIT  CBSEFVATICNS : 

a.  Wave  Period  - ftsccrd  the  tine  in  seconds  for  eleven  (11)  wave  "crests"  to  pass  so*  stationary  point.  Elsvsn  "crests"  will  includs  tan  caplets 
wavra  (crests  and  troughs) . Tha  first  (1)  "Chest"  selected  for  observation  is  recorded  as  tine  zero  and  the  eleventh  (11)  “creet"  will  be  the  stop  or  cut 
tint.  Record  this  tine  in  seconds  in  the  spaces  provided. 

b.  Brewing  Wave  Height  - This  observation  is  based  solely  an  the  ju&pnent  of  the  cbeerver.  Natural  or  marinade  features  cn  the  shoreline  or  in  the 
surf  zone  whose  dimensions  are  known  ray  aid  in  judging  the  height  of  a wave.  Otherwise  the  observer's  best  estimate  will  be  sufficient.  Record  the  breaker 
height  to  the  nearest  tenth  of  a foot. 

c.  Breaker  Angle  - To  ctetermine  the  direction  from  which  the  waves  are  approaching  the  bead!  use  the  protractor  an  this  reverse  side  of  the  data  form. 
The  0-180"  line  should  be  oriented  alcng  the  shoreline;  use  the  protractor  to  site  the  direction  from  which  the  waves  are  approaching  when  they  are  first 
breaking. 

d.  Type  of  Breaking  Wave: 

Spilling  - Spilling  occurs  wnen  the  wave  crests  becomes  unstable  at  the  top  and  the  crest  flews  down  the  front  face  of  the  we'*  producing  an  irregu- 
lar, foamy  water  surface,  (see  figure  1) 

Plinging  - Plunging  occurs  when  the  wave  crest  curls  over  the  front  face  of  the  wave  and  falls  into  the  base  of  the  wave  proAmng  a high  splash  and 
much  foon  (figure  2) 

Surging  - Surging  occurs  when  the  wave  crests  remains  unoroken  while  the  base  of  the  front  face  of  the  wave  adcanoee  up  the  beach [anr  figure  3T  - 
Spill/Plunge  - A certoinaticn  of  boch  spilling  and  plunging  oacurnnq  simultaneously. 

WIND  OBSERVATIONS : 

a.  Wind  Speed  - A wind  meter  is  provided  to  each  observer  and  it  is  recamended  that  the  instructions  provided  with  the  meter  be  followed  to  obtain 
wind  speed  measurements. 

b.  Wind  Direction  - After  the  approximate  orientation  of  the  beach  with  respect  to  north  has  been  defined  the  observer  can  determine  the  direction 
"from  which"  the  wind  is  coning. 

PORES HOPE  SLOPE: 

For  measurement  of  the  foreshore  one  rust  use  either  the  clipboard/ inclincreter  or  the  Abney  hand  level.  Cbservaticns  should  be  made  as  close  to  mid-*wash 
as  possible.  Using  the  clipboard/ inclinaneter  place  it  an  the  appropriate  edge  and  record  the  angle  where  the  ball  acmes  t*.  rest.  Using  the  Abney  .rand 
level  place  it  cn  a straight  edge  and  level  the  bifcble;  record  the  indicated  angle. 

WXDIH  CT  Sira  ZCNE: 

This  observation  is  based  solely  on  the  judgment  of  the  observer.  Estimate  in  feet  the  distance  from  the  shoreline  to  the  line  of  the  most  seaaard  breakers 
(not  to  be  confused  with  white  cape) . 

LCNCSHCBE  CURFOJT: 

a.  Dye  Distance  - Dye  packet  should  be  injected  just  shoreward  of  the  breakers,  if  possible.  Driftwood  or  any  other  floating  abject  shcxild  be  used  if 
dye  is  not  available.  Estimate  the  distance  from  the  shoreline  to  point  of  injection  and  record  this  distance  in  feet. 

b.  Current  Speed  - Mark  the  beach  in  line  with  the  injected  dye  and  make  a second  mark  to  indicate  the  dye  movement  after  cn*  minute  hae  lapsed.  Pace 
the  distance  between  these  marks  and  record  this  distance  in  feet. 

c.  Current  Direction  - When  locking  seacard,  if  the  dye  has  nwed  to  the  left  record  -1,  to  the  right  record  +1,  and  no  longshore  movement  record  0. 
RIP  anraNIS: 

Rip  currents  are  defined  as  seaward  moving  dvamels  of  water  which  return  the  water  that  has  been  piled  up  along  the  shore  by  incoming  waves.  Rip  currents 
are  fed  by  feeder  currents , water  moving  alona  the  shore  1 see  f igure  4) . TWo  currents  join  and  extend  out  in  what  is  known  as  the  "neck" , where  the  water 
rushee  through  the  breaker  zone  in  a narrow  lane.  Beyond  the  breaxer  zone  the  current  spreads  cut  and  dissipates  in  what  is  called  the  "head".  If  sutfr 
np  currents  are  present  estimate  their  spacing  in  feet.  If  no  rips  are  present  record  0. 

BEACH  CUSPS: 

Cusps  are  samcircular  or  crescent  shaped  cutouts  in  the  beach  face  see  figure  5) . If  such  shapes  are  observed  record  the  distance  between  the  "herns"  of 
the  cueos  which  indicate  the  spacing.  Where  the  spacing  is  irregular  estimate  the  average  spacing.  If  no  cusps  are  present  record  3. 


FIGURE  1 


SPILLING  WAVE 


FIGURE  2. 
PLUNGING  WAVE 


. — - 


l vr~ 


SUtRlC  FROM 


FIGURE  3 


SURGING  WAVE 


FIGURE  4. 


. • \> 

FIGURE  5 BEACH  CUSPS 


- - - . . - 


51 


Abramowitz,  M.  and  Stegun,  !.  A.,  Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards,  1964. 

Balsillie,  J.  H.,  Surf  Observations  and  Longshore  Current 
Prediction,  Technical  Memorandum  No.  58,  U.  S.  Army, 
Coastal  Eng . Research  Center,  November  1975. 

Battjes,  J.  A.,  Computation  of  Set-up,  Longshore  Currents, 
Run-Up  and  Overtopping  Due  to  Wind-Generated  Waves,  I 9 

Bowen,  A.  J.,  "The  Generation  of  Longshore  Currents  on  a 

Plane  Beach",  Journal  of  Marine  Research,  v.  37,  p.  206 
15,  1969. 

Bretschne ider,  C.  L.,  Field  Investigations  of  Wave  Energy 
Loss  of  Shallow  Water  Ocean  Waves,  U.  S.  Army,  Beach 
Erosion  Board,  Technical  Memorandum  46,  1954. 

Collins,  J.  I . , Long  Shore  Currents  and  Wave  Statistics  in 
the  Surf  Zone,  Technical  Report  No.  TC-149-2,  Tetra 
Tech,  Inc.,  February  1972. 

Galvin,  C.  J.  Jr.,  and  Savage  R.  P.,  Longshore  Current  at 
Nags  Head,  North  Carolina,  U.  S.  Army  Coastal  Engineer 
ing  Research  Center  Bulletin  II  (1966),  p.  11-29. 

Galvin,  C.  J.,  and  Nelson,  R.  A.,  Compilation  of  Longshore 
Current  Data,  U.  S.  Army  Coastal  Engineering  Research 
Center,  Miscellaneous  Paper,  March  1967. 

Galvin,  C.  J.  Jr.,  "Longshore  Current  Velocity:  A Review 
of  Theory  and  Data",  Review  of  Geophysics,  v.5.  No.  3, 
p.  287-304,  August  1967. 

Ingle,  J.  C.  Jr.,  The  Movement  of  Beach  Sand,  Elsevier 
Publishing  Company,  1966. 

Inman,  D.  L.  and  Quinn,  W.  H.,  Currents  in  the  Surf  Zone, 
in  Proceedings  2nd  Conference  on  Coastal  Engineering, 
edited  by  J.  W.  Johnson,  Council  of  Wave  Research, 

San  Francisco,  California,  p.  24-36,  1952. 

Komar.  P.  D..  Beach  Processes  and  Sedimentation,  Prentice- 


Longuet-H i gg i ns , M.  S.,  "Longshore  Current  Generated  by 
Obliquely  Incident  Sea  Waves",  Parts  I and  2,  Journal 
of  Geophysical  Research,  v.  75,  No.  33,  November  1970. 

Munk,  W.  H.,  "The  Solitary  Wave  Theory  and  Its  Application 
to  Surf  Problems",  Annals  of  the  New  York  Academy  of 
Sciences,  v.  51,  p.  376-462,  1949. 

Miller,  I.  and  Freund,  J.  E.,  Probability  and  Statistics 
for  Eng i neers , Prentice-Hall,  Inc.,  1965. 

Nie,  N.  H.,  and  others.  Statistical  Package  for  the  Social 
Sciences,  2nd  ed.,  McGraw  Hill,  1970. 

Prandtl,  L.,  Essentials  of  Fluid  Dynamics,  Haffner,  New 
York,  1952. 

Sonu,  C.  J.,  Computer  Prediction  of  Nearshore  and  Surf-Zone 
Stat  i sties,  Tetra  Tech,  Inc.,  September  1975. 

Thornton,  E.  B.,  Longshore  Current  and  Sediment  Transport, 
Technical  Report  No.  5,  Prep,  under  Cont.  No.  WP  00889, 
Dept,  of  the  Interior,  Federal  Water  Pollution  Control 
Administration,  December  1969. 

Thornton,  E.  B.,  Kinematics  o*  Breaking  Waves,  Proceeding 
of  the  15th  Coastal  Engineering  Conference,  American 
Society  of  Civil  Engineers,  461-476,  July  1976. 


INITIAL  DISTRIBUTION  LIST 


No. 


Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

L i b ra ry  (Code  0142) 

Naval  Postgraduate  School 
Monterey,  California  93940 

Department  Chairman,  Code  68 
Department  of  Oceanography 
Naval  Postgraduate  School 
Monterey,  California  93940 

Dr.  Edward  B.  Thornton,  Code  68Tm 
Department  of  Oceanography 
Naval  Postgraduate  School 
Monterey,  California  93940 

CDR.  Carlos  Severino  Veitia  Garcia 
Embajada  de  Venezuela 
Ag regadu r i a Naval 
2409  California  St.,  N.W. 
Washington,  D.  C.  20008 

Department  of  Oceanog ra phy , Code  68 
Naval  Postgraduate  School 
Monterey,  California  93940 

Oceanographer  of  the  Navy 
Hoffman  Building  No.  2 
200  Stova I I Street 
Alexandria,  Virginia  22332 

Office  of  Naval  Research 
Code  4 I 0 
NORDA,  NSTL 

Bay  St.  Louis,  Mississippi  39520 

Dr.  Robert  E.  Stevenson 
Scientific  Liaison  Office,  0NR 
Scripps  Institution  of  Oceanography 
La  Jolla,  California  92037 

Dr.  Doug  las  L.  I nman 
Scripps  Institution  of  Oceanography 
University  of  California 
La  Jolla,  California  92037 


Copies 

2 

2 

I 

5 

I 

3 


55 


11.  Library,  Code  3330 

Naval  Oceanographic  Office 
Washington,  D.  C.  20373 

12.  S I 0 Library 

University  of  California,  San  Diego 

P.  0.  Box  2367 

La  Jolla,  California  92037 

13.  Department  of  Oceanography  Library 
University  of  Washington 
Seattle,  Washington  98105 

14.  Department  of  Oceanography  Library 
Oregon  State  University 
Corvallis,  Oregon  97331 

15.  Commanding  Officer 

Fleet  Numerical  Weather  Central 
Monterey,  California  93940 

16.  Commanding  Officer 

Naval  Environmental  Prediction 
Resea  rch  Fac i I i ty 
Monterey,  California  93940 

17.  D i recto  r 

Naval  Oceanography  and  Meteorology 
National  Space  Technology  Laboratories 
Bay  St.  Louis,  Mississippi  39520 

18.  N0RDA 

Bay,  St.  Louis,  Mississippi  39520 

19.  Dr.  William  S.  Gaither 

Dean,  College  of  Marine  Studies 
Robinson  Hall 
University  of  Delaware 
Newark,  Delaware  39520 

20.  Department  of  the  Navy 

Commander  Oceanographic  System  Pacific 
Box  1390 

FP0  San  Francisco  96610 

21.  Coastal  Studies  Institute 
Louisiana  State  University 
Baton  Rouge,  Louisiana  70803 

22.  Mr.  Thorndike  Saville,  Director 
Coastal  Engineering  Research  Center 
Department  of  the  Army 

Fort  Belvoir,  Virginia  22060 


56 


